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Experiments on the iron-pnictide superconductors appear to show some materials where the 
ground state is fully gapped, and others where low-energy excitations dominate, possibly indicative 
of gap nodes. Within the framework of a 5-orbital spin fluctuation theory for these systems, we 
discuss how changes in the doping, the electronic structure or interaction parameters can tune the 
system from a fully gapped to nodal sign-changing gap with s-wave (Ai g ) symmetry (a ). In partic- 
ular we focus on the role of the hole pocket at the (n, n) point of the unfolded Brillouin zone identified 
as crucial to the pairing by Kuroki et ali^, and show that its presence leads to additional nesting of 
hole and electron pockets which stabilizes the isotropic s state. The pocket's contribution to the 
pairing can be tuned by doping, surface effects, and by changes in interaction parameters, which 
we examine. Analytic expressions for orbital pairing vertices calculated within the random phase 
approximation fluctuation exchange approximation allow us to draw connections between aspects 
of electronic structure, interaction parameters, and the form of the superconducting gap. 
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I. INTRODUCTION 

In any new class of superconductors, the structure of the order parameter is an important clue to the nature of 
the pairing mechanism, but the determination of this structure is seldom immediate. In the cuprates, for example, 
many different types of experiments on a variety of samples had to be analyzed and compared before a consensus 
was achieved, and the experimental picture was not clarified until the effects of disorder were understood and clean 
samples were prepared. It is therefore not unexpected that the symmetry and momentum structure of the gap in 
the Fc-based superconductors are still controversial nearly two years after their discovery— Nevertheless, the range of 
behaviour seen in different materials is striking^ Here we pose, from the point of view of a weak coupling fluctuation 
exchange theory^^ the questions: is it possible that the superconducting state of the ferropnictides is intrinsically 
sensitive to aspects of the electronic structure which "tune" the pairing interaction? If so, which degrees of freedom 
are most important? 

Based on density functional theory (DFT)^r— quantum oscillations and angle- resolved photocmission experiments 
(ARPES)2r— , the Fermi surface of the Fe-pnictides is believed to consist of a few small hole and electron pockets, as 
shown in Fig. [TJ where we have also indicated the predominant Fe-orbital character of the various parts of the Fermi 
surface taken from the DFT calculations of Cao et al& for the LaFeAsO material. We will follow the convention in 
Ref. H and elsewhere and refer to the hole pockets around the (0,0) point as the a sheets and the electron pockets 
around the (it, 0) and (0, it) points of the unfolded (1-Fe) Brillouin zone as the (3 sheets. Early on, it was proposed by 
Mazin et al^ and Dong et al^ that the nested structure of this Fermi surface would lead to a peak in the magnetic 
susceptibility near (tt,0), and that this might drive a sign change in the superconducting order parameter between 
the a and /3 sheets^. Several experiments on the Fe-based superconductors are indeed consistent with a gap which 
is isotropic (independent of momentum on a given pocket), but possibly with overall sign change of this type. Angle- 
resolved photoemission spectroscopy (ARPES) experiments, while not sensitive to the sign of the gap, are the most 
direct measure of its magnitude, and have consistently provided evidence taken to support an isotropic gap structure 
in momentum space^^ The observation of a resonance in inelastic neutron scattering is strong evidence for a sign 
change of the superconducting gap^~— . 



n = 6.01 n = 5.95 




FIG. 1. Fermi sheets of the five-band model for n = 6.01 (left) and n = 5.95 (right) with colors indicating majority orbital 
character (red=d xz , green=d HZ , b\ue=d xy ). Note the 7 Fermi surface sheet is a hole pocket which appears for ~ 1% hole 
doping. 

On the other hand, many other experiments appear to support the existence of low-lying excitations below the 
apparent gap energy. For example, in both the LaFePO syste m 23 i 24 and in BaFc2Asi_ a; P a f^; a linear-T dependence 
of the low-T penetration depth AA(T) has been reported, and in Ba(Fei_ a ,Co a ,)2As2, AA varies close to T 2 over 
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most of the phase diagram^ 6 .; these power laws are in contrast to the activated temperature dependences expected 
for an isotropic gap. Similar power laws have been observed in NMR21~— , thermal conductivit y 2 ^ 33 " — , and Raman 
scattering 3 ^ One obvious way of interpreting these observations is to suppose that the superconducting gap has nodes 
on parts of the Fermi surface, allowing for the excitation of quasiparticles at arbitrarily low energies. However, one 
may also show that in an isotropic "sign-changing s-wave" (s^) superconductor, disorder can create subgap states^ 
under certain conditions, depending on the ratio of inter- to intraband impurity scattering. An impurity band at 
the Fermi level in an state will also lead to AA ~ T 2 . There is no known scenario for producing A A ~ T with 
impurity scattering in a gapped state, however. It is extremely important to establish whether low-energy excitations 
are intrinsic (nodal) or extrinsic (disorder- induced) , and under what circumstances fully developed gaps, as opposed 
to highly anisotropic gaps, possibly with nodes, should be expected. 

From the standpoint of fluctuation exchange theories of pairing based on realistic Fermi surfaces in these materials, 
the most likely states indeed appear to be preferentially of " s-wave" symmetry, with quasi-isotropic gaps on the hole 
pockets but highly anisotropic states on the electron pockets^^. All of these calculations indicate the proximity of 
other pairing channels, particularly one with d x 2_ y 2 symmetry, but transitions between an s-wave state and a d-wave 
state would give rise to thermodynamic anomalies which have not yet been convincingly observed. Attention has 
therefore focussed primarily on the possibility of s-wave {A\ g symmetry) states with "accidental" nodes, i.e. nodes 
whose existence is due to details of the pairing interaction rather than symmetry. When the leading instability was 
of s-wave (Ai s )type, 5-orbital calculations found highly anisotropic states for all values of parameter space explored^, 
in apparent contradiction to the existence of nearly isotropic states experimentally observed in some materials. 

What aspects of the physics of these materials are responsible for the nodes or near-nodes seen in these theories? 
Some observations on this question have already been made. Maicr et al^ pointed out that, within a model with 
intra- and inter-orbital interactions, nodes were driven by the intra-orbital Coulomb repulsion, the scattering between 
the two /3 sheets neglected in simpler 2-band approaches, and a tendency (observed for the parameters considered 
in that work which were consistent with local spin rotational invariance) of electrons in like orbitals to pair. Kuroki 
et alX made an important connection between the lattice structure, electronic structure, and pairing state of the 
Fe-based superconductors, observing that in DFT calculations the pnictide atom height above the Fe plane appeared 
to control the appearance of a third 7 Fermi surface sheet centered on the T point in the folded zone corresponding to 
a (7r, 7r) pocket in the unfolded zone. This new hole-type pocket, not considered in Ref. [H, stabilizes a more isotropic 
state. When the 7 pocket is present, intra-orbital q ~ (-7T, 0) and (0,7r) scattering of d xy pairs between the 7 and 
(3 pockets favor a nodeless state. 

Within a model with band interactions, Vildosola et a/4 3 - and Calderon et a/4i have also discussed the change in 
the electronic structure caused by the shift of the pnictogen. In particular, the latter authors have noted that a change 
in the angle a formed by the Fe-As bonds and the Fe-plane can modify the orbital content as well as the shape of the 
Fermi surface sheets. In a similar model, Chubukov et al£± deduced a phase diagram manifesting a transition between 
a nodal and fully gapped state with s symmetry as a function of a parameter controlling the relative importance of 
intraband repulsion, and Thomale et aL— reached similar conclusions within a 4-band model, exploring the stability 
of the nodeless state with respect to doping and other changes in electronic structure. Ikeda et al^ considered 
the doping dependence of spin fluctuations and electron correlations within the renormalized fluctuation exchange 
(FLEX) approximation using LDA dispersions; they also observe anisotropic behaviour of the gap function around 
the M point for the electron doped material, and fully-gapped behaviour in the hole-doped material. In a subsequent 
work, they explored the doping dependence of the pairing state within FLEX42. Wang et alM- have also discussed the 
important role played by the 7 Fermi surface and emphasized the role of the orbital matrix elements in determining the 
momentum structure of the gap. From their functional renormalization group calculations, they find that the degree 
of gap anisotropy depends upon the relative weight of the orbital matrix elements. They find that by tuning some of 
the tight-binding parameters to alter the relative orbital weights, one can change the anisotropy of the gap. Thomale 
et al^ have also recently reported functional renormalization group results and argue that a nodal superconducting 
gap can appear on the j3 electron Fermi surfaces when the 7 Fermi surface is absent. 

Many authors have considered a Hamiltonian with all possible two-body on-site interactions between electrons in 
Fe orbitals as a good starting point for a microscopic description of this system, 

H = Hq + n i£t n in + U' X! n ^ n ^' C 1 ) 

i,l i,l'<i 
i,£'<£ o,a' 

where nu = l»,*t + n ni- The interaction parameters U, U', J, and J' correspond to the notation of Kuroki et alA, 
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and are related to those used by Graser et alJ± by U = U, J = J/2, U' = V + J/4, and J' — J' . The kinetic energy 
H is described by a tight-binding model 5 - spanned by the 5 Fe d orbitals, which, depending on the tight binding 
parameters and the filling n give rise to the Fermi surfaces shown in Fig. [TJ Here in the 1 Fe per unit cell Brillouin 
zone that we will use, there are the usual a\ and a 2 hole Fermi surfaces around the T point, the two j3\ and Fermi 
surfaces around (tt, 0) and (0,7r) and for the hole doped case n = 5.95 shown in Fig. [T] there is an extra hole FS 7 
around the (tt, tt) point. We use a notation in which a given orbital index I 6 (1, 2, . . . , 5) denotes the Fe-orbitals 
(d xz ,d yz ,d X y,d x 2_ y 2,d 3z 2_ r 2). An important role will be played by the orbital matrix elements a£(k) = (^fk) which 
relate the orbital and band states. The dominant orbital weights af,(k)| 2 on the Fermi surfaces are illustrated in 
Fig.ffl 

In Eq. ((2J, we have separated the intra- and inter-orbital Coulomb repulsion, as well as the Hund's rule exchange 
J and "pair hopping" term J' for generality, but note that if they are generated from a single two-body term with 
spin rotational invariance they are related by U' = U — 2 J and J' = J. Below we also consider the case where spin 
rotational invariance (SRI) is explicitly broken by other interactions in the crystal, such that J' and U' can take on 
independent values. 

The analysis of individual orbital contributions to pair scattering can be a powerful tool to understand the influence 
of electronic structure on pairing. The basic picture which emerges from the spin fluctuation theories^ 5 - is straight- 
forward. The intra-orbital scattering of d xz and d yz pairs between the a and (3 Fermi surfaces by (tt,0) and (0,7r) 
spin fluctuations leads naturally to a gap which changes sign between the d xz /d yz parts of the a Fermi surfaces and 
the d xz /d yz parts of the j3\ and (3% electron pockets. However, intra-orbital scattering between the d xy portions of 
the /3i and P2 Fermi surfaces suppresses the gap in the d xy regions of the (3 Fermi surface and competes with the 
formation of an isotropic s gap. In addition, this anisotropy can also be driven by the need to reduce the effect of 
the on-site Coulomb interaction ! 42 ' 45 There can also be inter-orbital pair scattering, such that for example a d xz pair 
on ai scatters to a d xy pair on j3\ . These inter-orbital scattering processes which depend upon U' and J' are weaker 
than the intra-orbital processes for spin rotationally invariant (SRI) parameters where U' + J' = U — J and J' = J. 
However, as noted by Zhang et al£^, the interaction parameters in the solid will be renormalized and one may have 
non spin rotationally invariant (NSRI) interaction parameters leading to enhanced inter-orbital pairings. 

Here we explore these issues using results obtained from a 5-orbital RPA calculation to further examine the effects 
of the electronic structure, the doping and the interaction parameters on the structure of the gap. Analytical results 
for orbital components of the pairing vertices have allowed us to understand the factors which influence the structure 
of the gap, and in particular the transition from nodal to nodeless behavior. In Section [II] we introduce the effective 
pair scattering vertex, its fluctuation-exchange RPA form and the pairing strength functional that determines the 
momentum dependence of the gap function^ Results for the gap function g(k) obtained for a typical set of spin 
rotationally invariant (SRI) interaction parameters are discussed in Section IIII1 The strongest pairings occur in the 
A\ g channel and the discussion in this section focuses on the role of the doping and the presence of the 7 pocket in 
determining whether the gap is nodeless. The nodeless character is found to be related to the additional resonant 
pairing provided by the 7 sheet, and is associated with the strong tendency for pairing between like orbitals when the 
7 pocket is of d xy character. The origins of the tendency to pair in like orbitals are then explored in Sec. IIIII Non 
spin rotationally invariant (NSRI) interaction parameters are discussed in Section IIV1 where we examine the effect of 
the Hund's rule exchange J and the pair hopping J' on the pairing interaction and the nodal structure of the gap. 
Section [V] contains a further discussion of the effect of the orbital character of the 7 hole pocket as well as the possible 
role of surface effects in causing nodeless behavior. Our conclusions arc contained in Section [VIII and in the appendix 
we discuss some approximate forms for the pairing vertex. 



II. SPIN FLUCTUATION PAIRING 

As in Ref. 0, we analyze the effective pair scattering vertex T(k, k') in the singlet channel, 



r y (k,k') = Re 



J2 a£'*(k)o£'*(-k) x r, lW4 (k,k',a; = Q)a% (k')a^ (-k') 



where the momenta k and k' are restricted to the various Fermi surface sheets with k e Cj and k' 6 Cj. The orbital 
vertex functions F^ 2 ^ 4 represent the particle-particle scattering of electrons in orbitals £±,£4 into £2^3 (see Fig. [5]) 
and in the fluctuation exchange formulatio n 52 ! 53 are given by 



lu s x f PA (k k',^ + \u s - it/ c xr A (k - k', W )l> + il> 



(2) 
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where each of the quantities U s , U c , xii °t c - represent matrices in orbital space as specified in the appendix. Note 
that the X^ PA term describes the spin-fluctuation contribution and the X^ PA term the orbital (charge)-fluctuation 
contribution. 





Intra Inter Mixed 

FIG. 2. Top: pairing vertex re 1 e 2 i 3 e i defined in terms of orbital states £< of incoming and outgoing electrons. Bottom: 
representative examples of classes of orbital vertices referred to in the text: intra-, inter- and mixed orbital vertices. 

For the parameter regions we will discuss, the dominant contribution to the pairing comes from the 5 = 1 particle- 
hole exchange given by the first term in Eq. The forms of the interaction matrices U s and U c are given in the 
appendix. As illustrated in Fig. [2] there are intra-orbital, inter-orbital, and mixed-orbital pair scattering processes. 
The contributions of each to the total pair scattering vertex in Eq. are quite different. In particular, as 
discussed below, the orbital matrix elements for k and — k states on the Fermi surface favor pairs which are formed 
from electrons in the same orbital state. We therefore find that, in spite of the fact that the mixed-orbital scattering 
can be significant, its contribution to the pairing interaction is negligible. 

If one writes the superconducting gap A(k) as Ag(k), with g(k) a dimcnsionlcss function describing the momentum 
dependence of the gaps on the Fermi surfaces, then g(k) is determined as the stationary solution of the dimcnsionlcss 
pairing strength functional* 



%(k)] 



E, 



dku 



ij Jd : J7IkT JCj up(k') 



^ 5 (k)r i3 (k,k')ff(k') 



dk\\ 
i ?Ci v F (k) 



[3(k)] 2 



(3) 



with the largest coupling strength A. Here the momenta k and k' are restricted to the various Fermi surfaces keC, 
and k' G Cj and i^„(k) = |Vk-E,/(k)| is the Fermi velocity on a given Fermi surface. The eigenvalue A provides a 
dimensionless measure of the pairing strength. 
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III. SPIN ROTATIONAL INVARIANT CASE 

Below we discuss the essential physics of the gapped-nodal transition within the constrained interaction parameter 
subspace where spin rotational invariance is assumed. It is worthwhile recalling the main points of the argument for 
an isotropic s ± state: 1) that a repulsive effective interaction peaked near (tt, 0) would drive such a state provided 2) 
this interaction did not vary significantly over the small Fermi surface pockets*^ Our main points are as follows: 

• The largest pair scattering processes tending to stabilize an isotropic state are the mira-orbital scattering 
pairing vertices r aaaa (c.f. Fig. [2]), which are peaked near (ir,0). 

• The intra-orbital processes primarily affect the gap on the sections of the Fermi surface with the corresponding 
orbital character. The relative signs between various orbital sections are determined by subdominant intra- and 
inter-orbital processes. 

• The isotropic state can be frustrated by intraband Coulomb scattering and by pair scattering processes 
between the two electron sheets with q ~ (77, 7r), both of which favor nodes42 

• The 7 pocket of d xy character can overcome this frustration and stabilize the nodeless stated 

• If the Hund's rule coupling J is weak, the processes induced by the 7 pocket are not sufficient to eliminate the 
nodes. The Hund's rule exchange is necessary to overcome an attractive (7r, 0) interaction between d xz (d yz ) 
pairs and d xy pairs and drive a strong intra-orbital d xz and d yz repulsion. 

To illustrate these points, we begin by plotting gap functions <?(k) for two typical sets of SRI interaction parameters 
U = 1.3, J = 0.0,0.2 and two different fillings n = 6.01 and 5.95, as shown in Fig. [3] In the electron-doped case 
(?i = 6.01) the 7 Fermi surface around (7r,7r) is absent and one finds an anisotropic gap with nodes on the f3\ and /?2 
Fermi surfaces for zero and for finite J £ The doping dependence at fixed J of the gap g(k) versus angle around the Pi 
Fermi surface is shown in Fig. SJ The nodes arise in part because of the f3i — (3 2 pair scattering and other aspects of 
the interactions which frustrate the isotropic s ± stated In particular, the sign change of the gaps on the end regions 
of the (3 Fermi surfaces reduces some of the frustration which arises from the j3\ — P2 scattering^ In addition, this 
sign change acts to suppress the effect of the intra-band Coulomb repulsion. 

For the hole doped case (n = 5.95) and finite J, the gap function, as seen in Figs. [3J and HI is anisotropic but 
nodeless. As discussed in Ref. [l], the scattering between the 7 (ir, n) Fermi surface and the /?i and /?2 Fermi surfaces 
stabilizes the nodeless gap. While there is still frustration associated with the (3i — $2 pair scattering, the additional 
fti — 7 and 02 ~ 7 pair scattering processes overcome it. The presence of the pocket also increases the overall pairing 
strength, which can be seen by comparing the pairing eigenvalues in the Figure; this would correspond to an increase 
in the critical temperature. To see how these effects arise in more detail, in Fig. [5] we have plotted several of the orbital 
pairing vertices r^ 2 ^ 3 f 4 (q) for q along high symmetry directions in the Brillouin zone for both n — 5.95 and 6.01. 
The peak in T2222(q) near the X point arises from q = (ir, 0) scattering processes in which d yz pairs are scattered 
between the a and (3\ Fermi surfaces. In the presence of the 7 sheet, both 1^3333 (q) intra-orbital d xy pair scattering 
between fi\^ and 7 as well as Fi33i(q) and T2332(q) inter-orbital d xz ,d yz pair to d xy pair scattering between /3i t 2 
and 7 are also possible. The pair scattering from the /3 sheets to the 7 hole pocket, represented by T3333, provides 
a strong resonant stabilization of the nodeless state; we have verified that the nodal behaviour is recovered when it 
is neglected. As discussed previously, although mixed-orbital vertices such as T2233 are present, their contribution 
to Tij(k, k') is suppressed by the orbital matrix elements in Eq. ([2]), and the pairing is dominated by intra- and 
inter-orbital scattering. 

To further illustrate this, in Fig. [S] we have plotted the pairing interaction rjj(k, k') which determines the gap via 
Eq. ([3]). In these plots, one member of a (k', — k') pair is located at the wave vector denoted by a black circle in the 
figure. The plot shows the strength of the pairing interaction (k, k') associated with the scattering of this (k', — k') 
pair to a (k, — k) pair on the various Fermi surfaces. One sees that if the initial pair is located in a region which 
has predominantly d^-orbital character (top), the strongest scattering is to a pair (k, — k) on other d xz regions^ 
Similarly, there are strong d yz intra-orbital scattering processes which are obtained by rotating the figures by n/2. 
While inter-orbital scattering processes are also present, they are weaker for the parameters we have considered here, 
as seen, e.g. from the plot of T2332 in Fig. [5] In addition to the d yz intra-orbital scattering, there is strong d xy 
intra-orbital scattering between /?2 and 7 as well as between /3i and 7. Note that in Fig. [5] the pairing strength is 
anisotropic along the Fermi sheets, violating one of the key assumptions underlying the argument for an isotropic 
state. 

We return to the role of the Hund's rule exchange J in stabilizing a nodeless state. The main effect, i.e. the lifting of 
the nodes as J is turned on, is illustrated in Fig.[7J There are two reasons this occurs. First, as discussed in appendix 
lAl the intra-orbital T's are controlled by U and J (see Eqs lA7l and lA8|) . For a fixed U, increasing J drives the system 



n = 6.01, J=0.0 



n = 6.01, J=0.2 



71 








1 


?i=.3236 


( ) 

mm 




- 






! I - 

i 



0.1 

0.05 



-0.05 
-0.1 



n = 5.95, J=0.0 



n = 5.95, J=0.2 



- 



I 


A,=.4237 


() 


_ 




— 

I 



0.1 

0.05 



-0.05 
-0.1 






A,=.7521 

O 


••••••• 

I 


o 

I 



i 



0.1 

0.05 



-0.05 
-0.1 







71 







7C 



FIG. 3. The gap eigenfunctions g(k) for a spin rotationally invariant parameter set U = 1.3, J = 0.0,0.2, for dopings n = 6.01 
(top) and n = 5.95 (bottom). 



closer to the instability and thus enhances the leading peak in the RPA susceptibility. However, Eq. @ also contains 
the term \{U S + U c ), which involves only the bare interactions rather than the RPA-enhanced susceptibilities. This 
term contributes to the intra-band Coulomb interaction (which favors anistropic pairing on the electron sheets) and 
scales, in the dominant orbital channels, with U rather than J. Thus increasing J increases the importance of the 
X\ PA term in Eq. @ relative to the intraband Coulomb interaction favoring the isotropic state. 

Secondly, we find that inter-orbital d xy to d xz pair scattering, represented by 1^331, plays an important role in 
stabilizing the isotropic state when a Hund's rule exchange is present. If J = 0, Fig. [5] shows that the pair scattering 
has an attractive peak at (tt, 0) which tends to induce nodes on the /3 Fermi surface. If J > 0, on the other hand, this 
attractive peak changes sign, as discussed in Sec. IIVI and in the appendix, and a nodeless state is stabilized by the 




(]) around (3 1 sheet 

FIG. 4. The gap function g(k) on the /3i pocket for n — 5.95, J = 0.2 (red triangles) and n — 6.01, J = 0.2 (blue squares). Here 
the angle <j) is measured from the fc^-axis. 



dominant intra-orbital pairing vertices. 

Finally, we note that in Graser et al£, upon electron doping a d-wave solution overtakes the anisotropic s-wave 
one. Here we find that upon strong hole doping (~ 8%), the spin- fluctuation approximation also leads to a d-wave 
solution. 



IV. BROKEN SPIN ROTATIONAL INVARIANCE 

In leading order, the strength of the inter-orbital pair scattering is determined by U' and J'. For SRI parameters 
U' = U — J — J' the intra-orbital pair scattering tends to dominate the pairing interaction as we have seen in the 
previous section. However, as noted by Zhang et ali^i, the actual interaction parameters appropriate for the Fe- 
superconducting materials need not be SRI. In this case, inter-orbital pair scattering may play a more important role 
in determining the momentum dependence of the gap function <?(k). In addition, using non-SRI parameters we can 
separately explore the roles of the exchange coupling J and the pair hopping J' interactions on the structure of the 
pairing interactions and the gap. Here, for a filling n = 5.95, we hold U and U' fixed and examine the roles of J and 
J' on the nodal-gapped transition on the /3i pockets near (tt, 0). The discussion below applies also to the 02 pocket 
near (0,7r) with orbital states rotated by 90°. 

In Fig. IH1 we plot both the various relevant orbital pairing vertices along high symmetry directions in momentum 
space, as well as the leading gap function <?(k) on the /3i sheet. Beginning with the nodal case J = J' = 0, we can see 
that increasing J only weakly affects the momentum space structure of the gap. However, note that the ^(k) plotted 
in Fig. [9] is normalized and should be scaled by the corresponding eigenvalue A to obtain the true gap amplitude. 
Turning on J indeed increases the pairing eigenvalue A by causing an overall increase in the pairing vertices. The 
largest scattering is provided by T2222 near (jt, 0), which is driven up by increasing J; this can be understood by 
observing that the RPA-rcnormalized susceptibility is enhanced by increasing J (see Eq. IA7|) . This clearly enhances 
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FIG. 5. Orbital pairing vertices along high symmetry directions in q space for U — 1.3, and J = 0.2 for n = 5.95 (bottom) and 
n = 6.01 (top), spin rotation invariance assumed. Solid (green) line: 1^2222 (intra); dotted (blue) ^332 (inter); dashed (red) 
T2233 (mixed); dashed (black) 1^3333; dash-dotted (magenta) Ti33i (inter). Note that the vertical scales in the two panels are 
different. 



the strength of the gap on the parts of the /?i Fermi surface characterized by a strong d yz orbital content. And, as 
discussed above, the repulsive nature of the interaction forces the gap on the d yz sections of the /3i sheet to have 
the opposite sign of the gap on the d yz sections of the a.\ and CV2 sheets. This also applies to the Tim vertex and 
corresponding d xz sections of the Fermi surface. The gap on the remaining portions of the Fermi surface, which 
are of d xy character, are left to be determined by the largest other pair scattering with a d xy component, which for 
these parameters is Fi33i. As discussed in appendix A, the Hund's rule coupling J contributes to an increase in the 
intra-orbital repulsive scattering which favors a nodeless gap, while also (see Eq. IA10[) leading to a more negative 
ri33i scattering, which favors a nodal gap. The net result, as seen in Fig. [5Jd), is that the nodes remain for J = 0.1 
and J' = 0. Alternatively, for the case in which J' =0.1 and J = 0, as shown in Fig. [9J3), the nodes are "lifted." In 
this case, J' is sufficient to overcome the negative contribution of U' 2 ^^- Here, has changed sign due to the 

contribution of Eq. IA9[ which is stronger than the contribution of U' (Eq. IA10[) due to its resonant structure. The 
change in sign now favors the same sign across the whole (5 sheet, which causes the nodes to lift. 
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FIG. 6. The total pair scattering vertex — k') for n — 5.95 with parameters U = 1.3 and J = 0.2 as a function of k with 

k' set to the point on the Fermi surface indicated in each panel by the black dot. 



V. ORBITAL CHARACTER OF HOLE POCKET 

We have discussed the appearance of the 7 hole pocket with doping in the context of doping by a rigid band shift 
of the 1111 Fermi surface of Cao et aZ.— . There appear to be various other scenarios in which electronic structure 
distortions might occur. For example, as noted by Kuroki et a/.—, variations in the As height, which are known to occur 
in the 1111 family, can tune the size of the 7 pocket and thereby the pairing itself, independent of doping. Another 
effect of the As height noted by Calderon et al£& is a switching of two bands near T such that the 7 pocket which 
occurs upon hole doping is of primarily c? 32 2_ r 2 character rather than d xy . Within our framework we can imitate this 
effect simply by adjusting the size of certain tight-binding coefficients associated with the Fe-As hopping in order to 
create such a d^ z 2_ r 2 pocket at 7, and ask what effect this has on the pairing. As also found by Calderon et al., the 
new 7 pocket which appears is the only Fermi surface sheet with d^ z 2_ r 2 character, so one may expect the pairing 
to be substantially altered relative to the situation with a 7 pocket of d xy character. In Fig. [TU1 we see that when 
the electronic structure is adjusted to create a d^ z 2_ r 2 pocket instead of a d xy , the pairing cigenfunction reverts to 
the nodal s-type found in the electron doped case, as expected since the additional 7-/3 condensation energy which 
stabilized the isotropic state has been lost. In addition, the pairing strength A is substantially reduced. 



VI. EFFECT OF SURFACE ON PAIR STATE 

As noted above, the presence of a hole pocket of mainly d xy character around (7r,7r) in the unfolded zone causes a 
nodeless state to be favored over a nodal one (c.f. Fig. [3]). As pointed out by Kuroki et alA, this can be accomplished by 
an increase of the pnictogen height /ip n . To provide some additional insight on the emergence of nodeless behavior, in 
particular in ARPES experiments^—, we have performed first-principles calculations on a slab of BaFe2As2 containing 
6 FeAs layers. 

The calculations were performed using the Quantum-ESPRESSO package^, which uses a plane wave basis. We used 
the Pcrdcw-Burke-ErnzcrhoS^ exchange-correlation functionals and ultrasoft pseudopotcntials. The use of ultrasoft 
pseudopotentials enabled us to utilize an energy cut-off of 40 Ry for the plane wave basis, while the density cut-off 
was taken to be 400 Ry. To determine the positions of the surface ions, two layers in the middle were kept fixed, and 
the rest of the atoms were allowed to relax. We find that the pnictogen height near the surface and the c-axis lattice 
constant contract, by about 13% and 5%, respectively. 

Fig- EH shows the band structure obtained for the BaFe2As2 slab. To determine the origin of the bands, we projected 
the band structure on the atomic wavefunctions of the ions. A particular is considered to belong to a certain ion if 
the projection onto the atom is larger than 50%. We have verified that the results below do not change appreciably 
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FIG. 7. The gap function g(k) on the /3i pocket for n — 5.95, J = (black diamonds) and n = 5.95, J = 0.2 (red triangles). 
Here the angle cj> is measured from the fc^-axis. 



if the projection threshold is varied. 

As can be seen from the figure, the presence of the surface causes the band just below the Fermi energy in the bulk 
to rise slightly, causing the appearance of an additional Fermi surface. Note that these results are in the folded zone; 
zone folding causes the pocket at (7r,7r) to appear at T, and it is this pocket that is seen due to the surface. For this 
pocket to cause a nodeless gap it is necessary that it is of d xy character (or d x 2_ y i in the folded zone), which we have 
confirmed (not shown). Thus, it is possible that due to the sensitive nature of the band structure, the presence of a 
surface can cause surface probes to detect a nodeless gap, even when the bulk gap has nodes. 



VII. CONCLUSIONS 



One of the striking features of the Fe superconductors is the sensitivity of the momentum space structure of the 
superconducting gap to relatively small changes in the electronic structure. Indeed, the electronic structure of these 
materials is quite special. First, they are semi-metallic with multiple, nearly compensated electron and hole Fermi 
surfaces. Secondly, multiple Fe orbitals lie near the Fermi energy. This means that relatively small changes in the 
doping or atomic structure can alter the nesting, the orbital composition of the band states on the Fermi surface, and 
even the number of k z =0 Fermi surfaces. The most prominent example of this phenomenon is the appearance of the 
so-called (vr, ir) hole pocket, studied first by Kuroki et ai, with hole doping or pnictogen height^ 

Within the framework of an RPA fluctuation exchange approximation, we can understand how these changes affect 
the structure of the superconducting gap. For the parameter regime that we have studied, the leading pairing state has 
A\ g symmetry. It is basically an s state, but the anisotropy of the gap, and particularly the presence or absence of 
gap nodes on the electron Fermi surface sheets is sensitive to the electronic structure. It is important to note that this 
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FIG. 8. The inter-orbital pair scattering vertex Fi33i(q) along high symmetry directions for n = 5.95 with parameters U = 1.3 
and J = 0.0,0.2. 



type of a unconventional nodal state is much more sensitive than for example, the d x 2_ y 2 pair state of the cuprates, 
where nodes owe their existence to symmetry; instead, in the Fe-pnictides the nodes appear to be "accidental" , i.e. 
determined by details of the pair interaction. 

As discussed in earlier works, the dominant pairing processes involve intra-orbital scattering. In this case, the 
intra-orbital (w, 0) and (0, it) scattering processes lead to a change of sign between the regions on the a and j3 Fermi 
surfaces where the d xz and d yz orbital weights are dominant. However, /3i to /?2 intra-orbital d xy scattering tends to 
frustrate the isotropic state. Furthermore, the effect of the intraband Coulomb interaction can be reduced in an 
anisotropic state. Thus, for light electron doping (n = 6.01), where there are only the a and j3 Fermi surfaces, we 
find that the gap has nodes on the /3 sheets. This anisotropy is further enhanced by the inter-orbital scattering of d xz 
(d yz ) pairs on a\ (a^) to d xy pairs on the /3 sheets. 

In the hole doped case (n = 5.95), or if the band structure is adjusted, a hole pocket appears around the point. 
For our band structure, this 7 pocket has d xy character. Intra-orbital d xy scattering between the j3 and 7 pockets 
favors a more isotropic gap, removing the nodes on the (3 Fermi surfaces. We note that if the band parameters were 
such that the orbital character of the 7 Fermi surface were, for example, d 3z 2_ r i instead of d xy , the gap nodes would 
return to the (3 Fermi surfaces. In addition, the pairing strength A is substantially reduced, reflecting the important 
role played by the orbital weights on the Fermi surfaced This illustrates the important role played by the orbital 
weight. 

We also pointed out that various effects in addition to the pnictogen height identified as a key factor by Kuroki et al. 
can influence the 7 pocket. In particular, the presence of a surface can create such a pocket in a nominally electron- 
doped system, stabilizing a nodeless state. This may explain why ARPES experiments have reported quasi-isotropic 
gaps in these systems. 

We have investigated in some detail the way in which the various single-site interaction parameters U, U', J and 
J' influence the different types of pair scattering processes. The RPA technique is crude, but it allows analytical 
insights into questions of this type. Wc have been therefore able to trace the effects of varying these parameters to the 
strength of the relevant orbital pair scattering vertices at appropriate nesting vectors, which in turn determines not 
only the overall magnitude of the pairing strength (i.e., T c ), but also the anisotropy of the gap on the Fermi surface. 
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FIG. 9. Left: orbital pairing vertices, and right: gap function on the /3i sheet, for U = 1.2, U 1 = 1.2 and n = 5.95. Cases 
shown are J = J'=0 (a, A = 0.10); J = 0.1, J' = (b, A = 0.21); J = 0, J' = 0.1 (c,A = 0.41). 



It is to be hoped that the influence of various changes in crystal structure, morphology, etc. on the pair state and 
transition temperature may now be better understood through this type of analysis. 
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Appendix A: The scattering vertices 



The basic scattering vertices for the multiorbital Hubbard model are shown in Fig. IA.11 Here, as noted in the 
introduction, we use a notation in which an orbital index I = (1, 2, . . . , 5) denotes the Fc-d orbitals (d xz , d yz , d xy , 
d x 2_ y 2, d 3z 2_ r 2). As seen, in lowest order, there are intra-orbital (a), inter-orbital (d) and mixed-orbital (b) and 
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n = 5.98, J=0.2 




FIG. 10. The pairing eigenfunction for the hole-doped (x=2%) compound where the pocket character has been adjusted to be 
of d 3z 2_ r 2 type. The interaction parameters have been chosen as U = 1.3, J = 0.2. 



(c) pair scattering processes. Some second order contributions are shown in Fig. IA.1I (e)-(g). The orbital matrix 
elements for k and — k states on the Fermi surface favor pairs which are formed from electrons in the same orbital 
state. Thus in spite of the fact that the mixed-orbital scattering can be significant, its contribution to the pairing 
interaction is negligible and the intra- and inter-orbital scattering processes give rise to the superconducting pairing^. 
The relative Z-orbital contribution to the Bloch state k on the v th Fermi surface is given by the square of the orbital 
matrix element |a^(k)| 2 . As shown in Fig. [TJ the I = 1 (d xz ) and I = 2 (d yz ) orbitals give the main contributions 
on the a Fermi surfaces. Likewise the / = 1 and I = 3 (d xy ) orbitals contribute to the fa Fermi surface, while the 
I = 2 and I = 3 orbitals contribute to the fa Fermi surface. For our tight-binding fit of the Cao et al. bandstructure^, 
the 7 Fermi surface (around M=(n, 7r)) has predominantly / = 3 (d xy ) weight. The orbital weight distribution favors 
(k, — k) pairs with similar orbitals so that both, intra-orbital (d xz ,d xz ) a -B- (d xz ,d xz )j3 2 , (d yz ,d yz ) a <H> (d yz ,d yz )/3 1 as 
well as inter-orbital (d XZl d xz ) a ■(-> (d xy ,d xy )i3 2 and {d yz ,d yz ) a <-> {d xln d xy )p 2 pair scattering processes contribute. In 
addition, when the 7 Fermi surface around (7r, 7r) is present there are important intra-orbital (d xy , d xy )j o {d xy , d xy )p 1 
and (d xy ,d xy ) 7 -o- (d xy ,d xy )i3 2 pair scatterings. From Fig. IA.1"! one sees that the first order intra-orbital scattering 
processes involve U and J while the inter-orbital processes depend upon U' and J'. 

In the RPA fluctuation exchange approximation, for the parameter range of interest, the dominant pairing interac- 
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FIG. 11. The DFT band structure calculated for a BaFe2As2 slab (gray points). The red points show the bulk contributions 
from the FeAs layers, while the black points denote the corresponding surface contributions. 



tion is given by Eq. © with 

xf PA (q) = X (q)[l-C/ s X°(q)] _1 (Al) 

and 

Xo RPA (q) = X°(q) [1 + C/ c X°(q)] ~* (A2) 
For the 5-orbital problem II s , U c and x° are 25 x 25 matrices in an (£±£2) basis with 




= £2 = £3 =£4 



w £ 1 =£ 3 ^£ 2 = £ i (A) 

£i=h^h = £3 

£x=h=h = U 
ii = h± £2 
h = £i + £3 



1 2J, £\ — £ 3 7^ £ 2 — £4 f \ a\ 

u tit2t,u - s ofr> lx = i2 -l t 3 = £ 4 y M > 



and 

1 ^ aft (k)<ff* (k)<# (k + q)^* ( k + q ) 



x° 1 £ 2 £ 3 £ 4 (q) 



N ^ £„(k + q)-.E M (k) 



x[/(£„(k + q))-/(£„(k))] (A5) 

Here \i and v are summed over the band indices and / is the usual Fermi function. We have typically taken the 
temperature T = 0.02 eV. The orbital indexing convention for the susceptibility is illustrated in figure [A~2l 

As seen in Fig.[IJ three orbitals 1 (d xz ), 2 (d yz ), and 3 (d xy ) have significant weight on the Fermi surfaces. Therefore 
in this case, the interaction matrix U s reduces to the 9x9 matrix shown in Table HI Furthermore, as seen in Fig. IA.31 
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//ft lk'\ //ft „ lk'\ //ft /'/f't //ft _ l'k'\ 





the bare off-diagonal elements of the susceptibility involve single-particle propagators projected on different orbitals 
which makes them smaller than the diagonal elements. 

Keeping only the diagonal terms in xo, it is straightforward to evaluate the RPA susceptibility matrix. For example, 



\ RPA / o 

X1313 X1331 \ _ / Xl313 U 



X3113 X3131 / \ X3131 



(A6) 



1-^X?313 -J'X° 



1313 



-J'X°3131 1-^'X§131 
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11 22 33 12 21 13 31 23 32 



11 

22 
33 
12 
21 
13 
31 
23 
32 



U J J 
J U J 
J J u 



U' J' 
J' V 



U' J' 
J' U' 



U' J' 
J' U' 



TABLE I. Interaction matrix II s in the reduced 1 {d xz ), 2 (d yz ), 3 (d xy ) basis. 



<n> = 6.01 




X M 

FIG. A. 3. The non-interacting susceptibilities X < e 1 e 2 i 3 t4 f° r n ~ 6.01. 



Using the fact that X1313 = X313H and setting X1313 = X13: we have 



rpa = 1 ( 1 , 1 

X1313 2 Xl3 ll-(t/'lJ')x! 3 !-({/'- J')X°13 

RPA = 1 ( 1 1 

X1331 2 Xl3 ^ 1 _ {t)l + J, )x 3 l _ [t j,_ j l)x 3 
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In the same way, from the upper left 3x3 part of the interaction matrix one obtains for example 



_x?((i-tfx§) (i-UxD-Jxlxt) 



Aim — ^ \ AI ) 

where D is the determinant of 1 — U s x° and takes the form 

D = (l-UxD (1-^) (1-^) 

-^ 2 (l"^3)x?X^2J 3 x ? x o x o (A8) 

with X ° = Xaaaa- ^ is clear from these expressions that, for repulsive interactions and within the current approx- 
imation, while the intra-orbital spin susccptibities depend on U and J, the inter-orbital susceptibilities depend on 
U' ± J'. 

Although we have only kept the diagonal terms in the bare susceptibility, due to the structure of the interaction 
matrices there are off-diagonal terms in the RPA enhanced susceptibility. From these components, we focus on the 
intra- and inter-orbital pair scattering processes, as the mixed-orbital processes are suppressed by the external matrix 
elements. 

Within this diagonal bare susceptibility approximation, the inter-orbital pair scattering elements of T are simple 
because they involve only a 2x2 interaction matrix. For example, the pairing strength for the spin- fluctuation 
scattering of a d xz pair to a d xy pair is 



1 1331 - ~ A 



(U' + J') 2 X° 13 {U 1 ~ J') 2 X?3 



l-(U>+ J')x? 3 1-{U'- J')X°13 



(A9) 



There is a corresponding T2332 element which describes the inter-orbital scattering of d yz pairs from the a Fermi 
surfaces to the d xy regions of the electron Fermi surface. 

Since X13 is associated with the nesting of the d xz and d xy parts of ct\ and /?i respectively, it peaks for wave 
vectors near (it, 0). For reasonable values of the pair hopping J', there can be a significant inter-orbital repulsive 
pair scattering from the d xz regions of the hole Fermi surfaces a to the d xy regions of the /? Fermi surfaces. As J' 
goes to zero, the diagonal approximation for this inter-orbital d xz to d xy pair scattering vanishes linearly with J'. If 
J' becomes small so that the neglected off-diagonal component of the susceptibility matrix X1331 are large compared 
with |j7Xi3i3j then there will be off-diagonal corrections to Fx3i3. The leading correction is given by 

r 133 i * U'(U' + 2J) X ° 1331 (A10) 

which gives the attractive interaction at (tt, 0) for J = J' = 0, as shown in Fig. [5] In lowest order the (C^') 2 Xi33i 
contribution comes from the second order processes as seen in Fig. IA.ll As seen in Fig. IA.31 X1331 is negative near 
the X point where q ~ (ir, 0). 

The intra-orbital pair scattering involves a 3 x 3 matrix and is in general more complicated. However, in the diagonal 
susceptibility approximation, the intra-orbital T matrix will have the same determinant in the denominator as given 
by Eq. (|A8[) . From this, one can see how U and J enter in determining the intra-orbital pair scattering. In particular, 
the exchange J couples to the various bare susceptibilities such that, for example, the peak in X2222 nea r X is reflected 
in Tim and T3333 as well as T2222 when J is finite. To summarize, in the diagonal bare susceptibility approximation, 
the leading behavior of the intra-orbital pair scattering is controlled by U and J. Increasing these interactions increase 
the repulsive scattering at q ~ (ir, 0) and (0,7r). The leading behavior of the inter-orbital scattering is controlled by 
U' and J' and in the diagonal bare susceptibility approximation vanishes as J' goes to zero. When J' is small, the 
dominant contribution to the inter-orbital d xz to d xy and d yz to d xy scattering varies as U'(U' + 2J)x°33i- 



Appendix B: The orbital gauge 



In the following section we want to explicitly show that the choice of a gauge for the orbital basis does not affect 
the physical quantities, although one has to carefully take the orbital gauge into account when comparing orbital 
dependent interaction parameters, as e.g. the pair hopping term J' . First we define a gauge transformation of the 
orbital creation and annihilation operators as 

r -d — r-„ p^ 1 — c 1 P -i< ^ 
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where cj fcr creates a particle with spin a in orbital I at site i. Diagonalizing the non-interacting part of the Fourier- 
transformed Hamiltonian the orbital phases can be absorbed into the matrix elements in the form 

if 

where \jr vcr {k) creates now a particle with spin a and momentum k in band v and the matrix element a„(k) = e l ^ e a^(k). 
For the bare multiorbital susceptibility as defined in Eq. IA5I we thus find the following transformation 

while the spin susceptibility as an observable is gauge invariant 

Xs{q) = 2 $UW?) = Xs(q)- 

If we now proceed to the interaction Hamiltonian Eq. [5] we note that the pair hopping term is not gauge invariant 

f V r f r f r;nr- e ,*.- T V P 24 U<^+<M r t t r ..,. r . e ,. 

and we can define the gauge transformation for J' as 

e 2i(-<f> e +<f> e ,) J) 

Using this relation allows us to write the gauge transformed interaction matrices as 

fj s / c —Tj s / C f ,i(<j>i 1 -4>e 2 -<t>e 3 +<l>e i ) 

It is straightforward to show that this relation also yields the correct gauge transformation for the RPA enhanced 
multiorbital susceptibility 

xr A (q) = x°(q) + x? PA (q)u s x (q) 

with 

(x? PA )w 3 iAq^) = rtwAt!.^ 1 -*"^ 1 

under this transformation, the orbital dependent pairing vertex (Eq. [5]) transforms as 

f^ lW4 (fc,fc , ,a;)=r^ W4 (A; ) fc / ,a;)e- i ^-^-^3+^4) 

and we see that the effective pairing vertex entering the linearized gap equation is gauge invariant. 

In the above, we have presented the orbital pairing vertices (in Figs. [SJ [8] and [9]) and orbital susceptibilities (in Fig. 
IA.3|) . Both of these quantities depend on the choice of of orbital gauge. We have chosen to present them in the basis 
where (1) the d xz and d yz orbitals are aligned along the Fe-Fe directions, and (2) none of the orbitals have a purely 
imaginary phase with respect to any other. In the notation above, (j>t = 0, for all I In this basis, a\,(— k) = a l *(h), so 
that the sign of the intra- and inter-orbital pairing vertices accurately reflects their contribution to the gauge invariant 
pairing vertex r ?J (k, k'). 
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